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摘 要 以 河南 省 新 乡 市 娘娘 庙 污 灌区 一 块 约 328 m? 的 田 块 为 研究 区 ， 采 用 均值 、 变 异 系数 、 半 方差 函数 及 
分 维 数 等 相 结合 的 方法 ， 分 析 一 次 污水 灌溉 前 后 10 d 内 东西 向 剖面 A、B、C 和 D( 南 北向 每 相距 10.5 m 设 
1 个 剖面 ， 由 北 至 南 依次 编号 为 A、B、C 和 D) 的 土壤 有 效 锅 的 垂 向 空间 变异 及 其 分 布 特征 ,为 田 块 土壤 中 有 
效 态 锅 分 布 与 土壤 植物 修复 相 结 合 的 研究 思路 提供 基础 性 探讨 。 结 果 表 明 : 灌水 后 ， 土 壤 有 效 锅 含量 在 剂 面 
A、B、C 和 D 的 均值 分 别 减 少 0.06 mg.kg 1、0.11 mg:kg 1!、0.14 mg:kg ! 和 0.23 mg.kg-!， 呈现 中 等 离散 程度 ; 
在 0.393~1.308m 变 程 范围 内 ， 块 基 比 CW(Co+C)<25%, 土壤 有 效 锅 剖 面 总 体 表 现 出 由 随机 性 引发 的 变 蜡 较 小 , 空间 
相关 性 较 强 ; 剖面 A、C 和 DD 土壤 有 效 态 锅 的 分 维 数 降低 ， 空 间 变 异 程度 减弱 。Kriging 插值 结果 表明 该 区 有 效 包 
呈 条 带 状 分 布 格局 ， 且 在 0~15 cm 土 层 呈 现 累积 现象 ; 沿 灌水 方向 ,土壤 有 效 锅 含 量 剖 面 等 值 线 分 布 由 密集 到 稀 
朴 。 相 关 分 析 表 明 ,土壤 有 效 锅 含量 与 土壤 有 机 质 含量 呈正 相关 ,与 土壤 含水 率 呈 中 等 程度 相关 , 与 pH 呈现 强 负 
相关 ， 表 现 为 土壤 有 机 质 含量 越 高 土壤 有 效 锅 含量 越 大 , pH 越 高 土壤 有 效 锅 含量 越 小 ; 污水 灌溉 对 土壤 有 效 锅 、 土 
壤 含 水 率 、 土 壤 pH 和 土壤 有 机 质 的 影响 显著 ， 同 时 会 削弱 有 效 锅 与 各 影响 因子 间 的 相关 关系 。 伴 随 着 污水 灌溉 动 
态 的 有 效 锅 分 布 和 变异 特征 可 以 为 应 用 植物 吸附 进行 土壤 重金 属 污染 修 复方 面 的 研究 提供 基础 数据 。 
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Abstract There is limited development of agricultural economy in the north area of China that has been irrigated with 
sewage water for a long time, because of the uneven distribution of water resources and different levels of soil pollution. The 
main contaminants of sewage irrigation are heavy metals, which contamination area has been expanding over the years. A 
number of large-scale or multi-scale studies have been conducted on heavy metal pollution in surface soils in irrigated areas. 
However, studies have less focused on the vertical distribution of heavy metal, but mainly concentrated on the characteristics 
of migration and leaching of heavy metals in the soil. It is therefore necessary to focus on the vertical cross-sectional view of 
the spatial variability of soil heavy metals in sewage irrigated areas on a small scale. In this study, we investigated the spatial 
variability and distribution of available cadmium in the soil after once sewage irrigation. Using classical statistics and 
geo-statistical techniques, we measured the spatial variability and distribution in a 328 m” field in the Niangniangmiao sewage 
irrigation area in Xinxiang, Henan Province. The classical statistical analysis suggested that the mean available cadmium in 


Soil sections A, B, C and D (the south-north sections with 10.5 m away from each other and across the field from north to south) 
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decreased, respectively, by 0.06 mg-kg ', 0.11 mgkg ', 0.14 mgkg ! and 0.23 mgkg ! after sewage irrigation; suggesting a 
moderate degree of dispersion. The dispersion in Soil section B increased by 4.32%, while it decreased in other sections by 
9.39% (A), 6.12% (C) and 10.98% (D), all of which were within the 28.319%-45.16% range after sewage irrigation. 
Geo-statistical analysis showed that the overall variation in available cadmium in the soil sections was random，Co(Co+C) < 
25%. Soil available cadmium had a strong spatial correlation within the 0.393-1.308 m range. The optimal empirical 
semi-variogarm models for simulating soil cadmium content for soil sections A, B, C and D were respectively spherical, 
spherical, Gaussian and exponential modles. Soil profile B had the highest (1.997) fractal dimension of available cadmium. 
The fractal dimension decreased after sewage irrigation suggesting a weakening spatial variability of soil available cadmium. 
Structural factors played a decisive role in the vertical spatial variability and distribution of available cadmium in the soil 
sections. Kriging interpolation suggested a banded distribution pattern. There was significant available cadmium accumulation 
in the 0—15 cm soil layer. The distribution contour of available cadmium in the Soil sections varied 位 om dense to sparse along 
irrigation direction. Correlation and variance analysis showed that available cadmium content was positively correlated with 
organic matter content, moderately correlated with soil moisture, strongly negatively correlated with pH. The effect of 
irrigation on each factor was significantly, moreover, irrigation could dampen the impact of the factors on available cadmium. 
The results on distribution and variation in soil available cadmium showed leaching and migration of soil available cadmium 
after sewage irrigation. This was more visible for the 0—15 cm top soil than the 30-40 cm sub-soil. The spatial variability of 
soil available cadmium was mainly influenced by soil structural factors. The study provided the primary investigation for 
phytoremediation of heavy metal in polluted soils. 


Keywords Sewage irrigation; Soil profile; Available cadmium; Spatial variability; Heavy metal 


土壤 特性 的 空间 差异 性 "使 重金 属 在 土壤 中 具 
有 一 定 的 空间 分 布 特征 。 自 20 世纪 60 年 代 由 
Matheron 系统 化 提出 了 “区 域 化 变量 ”的 概念 路，80 
年 代 McBratney 等 中 形成 地 质 统计 学 理论 ， 到 现 阶 
段 GIS 和 Kriging 插值 的 综合 应 用 ,土壤 特性 空间 变 
异 定量 化 的 研究 得 到 了 较为 深入 的 发 展 ， 并 成 功 应 
用 于 研究 土壤 水 分 特征 S44、 水 盐 动 态 中 5 、 土 壤 养 


析 有 效 锅 与 其 他 影响 因素 间 的 相关 关系 ,为 该 污 灌 
区 土壤 中 有 效 态 锅 分 布 与 土壤 植物 修复 相 结合 的 研 
究 思路 提供 基础 性 探讨 。 


1 研究 区 概况 与 研究 方法 


1.1 ”研究 区 概况 
河南 省 新 乡 市 地 处 河南 省 北部 ， 黄 河中 下 游 ， 


分 "I、 重 金属 ”及 其 他 元 素 的 空间 分 布 上 。 我 国 
华北 地 区 ， 因 水 资源 分 配 不 均 ,利用 污水 灌溉 农田 
历史 悠久， 由 污水 灌溉 带 来 的 不 同 程度 的 土壤 污染 
制约 着 污 灌区 农业 经 济 的 发 展 。 重 金属 作为 污 灌 主 
要 的 污染 物 之 一 ， 影 响 面积 在 不 断 扩 大 ， 尤 以 锅 、 示 
及 其 复合 污染 较为 突出 中 I。 土壤 重金 属 在 不 同 地 区 
表现 出 不 同 的 变异 特征 ， 在 一定 范围 内 具有 明显 的 
空间 相关 性 中 "1 在 同一 地 区 不 同 尺 度 下 ,部 分 重 
金属 的 空间 估 值 可 以 较 好 地 表现 变量 空间 变异 1。 
然而 , 现 有 研究 大 多 是 关于 大 尺度 或 多 尺度 灌区 的 
耕 层 土壤 重金 属 污染 现状 J， 也 有 少数 学 者 对 土 
塘 重 金属 的 垂 向 分 布 进行 研究 ,但 主要 是 关注 重金 
属 淋 失 迁移 的 特征 。 长 期 灌溉 下 ， 重 金属 元 素 表 现 
出 在 土壤 表层 累积 ， 土 壤 锅 的 淋 失 迁移 程度 大 于 土 
壤 铅 中 。 因 此 ， 对 于 灌区 小 尺度 纵向 剖面 土壤 重 
金属 的 空间 变异 还 需要 进一步 研究 。 本 研究 以 河南 
省 新 乡 市 长 期 进行 污水 灌溉 的 田 块 作为 研究 区 域 ， 
运用 传统 统计 学 、 地 质 统计 学 、Kriging 插值 、 分 形 
理论 和 相关 分 析 等 方法 ， 探 讨 了 污水 灌溉 前 后 垂 向 
剖面 土壤 有 效 饥 含 量 的 空间 变异 性 与 分 布 特征 ,分 


为 驯 北 政治 、 经 济 、 文 化 和 交通 中 心 之 一 ， 也 是 国 
内 较为 重要 的 电池 生产 地 。 研 究 田 块 位 于 新 乡 市 郊 
娘娘 庙 灌 区 , 该 灌区 长 期 大 面积 采用 冬小麦 - 夏 玉 
米 的 耕种 模式 ,经 水 渠 引 卫 河 水 灌溉 。 卫 河 通 常 系 
指 河南 省 新 乡 市 合 河 乡 至 山东 省 聊城 市 征 钓 湾 一 段 ， 
至 西向 东 流 经 娘娘 庙 灌 区 ,年 均 流量 约 为 2 mi.s 1。 
卫 河 因 长 年 接纳 引 黄 灌溉 退 水 及 城市 生产 生活 污水 ， 
逐渐 成 为 生态 失控 的 排污 河道 1, 但 娘娘 庙 灌区 并 
未 因 卫 河水 质 逐 年 变 差 而 改变 其 灌溉 水 源 ， 可 见 利 
用 污水 灌溉 依然 是 缓解 当地 农业 用 水 紧缺 的 主要 途 
径 之 一 。 研 究 区 内 土壤 主要 为 粉 砂 质 黏土 ， 灌 溉 水 
中 ,总 氮 合 量 超出 《农田 灌溉 水 质 标准 》(GB 5084 一 
2005) 旱 作 标 准 值 ， 重 铬 酸 盐 指数 (CODcr) 和 总 氮 含 量 
超出 《地 表 水 环境 质量 标准 》(GB 3838 一 2002)V 类 。 
1.2 ”研究 方法 与 数据 处 理 
1.2.1 样本 采集 

本 次 研究 的 田 块 为 一 块 长 42 m( 南 北向 )、 宽 7.2 m 
(东西 向 ) 的 长 方形 区 域 , 考虑 到 取 土 、 实 验 室 检测 等 
的 工作 量 ,在 南北 向 每 相距 10.5 m 设置 1 个 剖面 ， 
由 北 至 南 依次 编号 为 A、B、C 和 D, 首尾 剖面 距离 
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地 块 边界 5.25 m; 东西 向 每 隔 0.9 m 布设 1 个 采样 点 ， 
则 每 个 剖面 上 为 7 个 采样 点 ， 首 末 采 样 点 与 田 块 边 
界 间距 0.9 m; 依据 种 植 小 麦 之 前 的 土壤 本 底 值 以 及 
重金 属 离子 在 土壤 中 迁移 缓慢 , 且 分 布 变化 不 明显 
等 特性 ， 取 土 过 程 使 用 钻头 直径 $S cm、 设 计 采 样 深 
度 为 1 m 土 钻 ,进行 0~40 cm 深度 的 土 样 采集 ， 以 
10 cm 为 间隔 。 每 个 剖面 上 以 长 方形 网 格 状 分 布 28 
个 采样 点 ， 测 定 污水 灌溉 前 后 土壤 含水 率 、pH、 有 
机 质 含量 和 有 效 态 锅 4 项 参数 。 

灌水 时 间 为 2013 年 4 月 25 日 ,两 次 采样 的 时 
间 分 别 为 2013 年 4 月 15 日 和 2013 年 5 月 2 日 , 采 
样 间 隔 期 间 未 出 现 降水 ; 灌水 量 约 为 60 mm 次",， 灌 
水 以 大 水 漫灌 方式 (灌溉 水 入 渗 深度 40~60 cm)， 垂 
直 于 剖面 方向 由 南 至 北 进行 灌溉 。 

试验 开始 时 ， 以 “S” 状 在 整个 地 块 取 5 个 采样 位 
置 ,每 个 采样 位 置 取 3 个 采样 孔 ,采样 深度 为 80 cm， 
以 10 cm 为 间隔 , 各 层 土 样 为 3 个 采样 孔 的 混合 
样 ， 经 检测 取 不 同 深度 的 平均 值 ， 田 块 土壤 基本 理 
化 性 状 为 : 锅 1.02 mg:kg :、 有 效 锅 0.356 mg.kg 、 
pH 8.5、 有 机 质 质量 百分数 2.14%、 容 重 1.42 gm ， 
0~40 cm 均 为 粉 砂 黏 质 士 ; 卫 河水 样 : pH 7.38、 电 导 率 
2 913.33 hnS.cm CODcr 72.77 mg 总 氮 23.17mg:IL、 
总 钾 12.87 mg.LT、 锅 0.004 969 mg:L"'、 铅 0.001 402 mg 
和 铬 0.002 107 mg:.L-L。 
1.2.2 ”样品 检测 

土 样 采集 过 程 中 , 一 钻 土 分 为 两 部 分 : 一 部 分 
湿 土 装 入 铝 盒 , 于 105 烘箱 中 烘 12h 后 取出 称 重 ， 
测定 土壤 含水 率 ; 另 一 部 分 湿 土 装 入 密封 袋 ， 于 实 
验 室 风干 ,研磨 , 过 2 mm 尼龙 科 , 用 DTPA( 二 乙 三 


胺 五 乙酸 ) 提 取 剂 以 1 : 2 制 成 悬 浊 液 ， 置 于 塑料 管 

， 震 荡 2 h 后 过 滤 得 清 液 PY， 使 用 日 本 岛 津 
(AA-6300) 原 子 吸收 分 光 光 度 计 测定 土壤 重金 属 镶 
的 有 效 态 含量 , 采用 低温 外 热 重 铬 酸 钾 氧化 - 比 色 
法 测定 土壤 有 机 质 ， 由 上 海 精 科 (752N) 柴 外 可 见 分 
光 光 度 计 测定 标准 曲线 和 各 土 样 的 吸光 度 值 ; 用 超 
纯 水 处 理 干 土 1 : 5 制 成 悬 浊 液 ， 采 用 电位 法 测定 土 
壤 pH, 使 用 pH 计 上 海 精 科 (PHS-3E) 酸 度 计 进行 测 
定 。 所 有 参数 指标 的 测定 方法 均 参 照 《 土 壤 农 业 化 
学 分 析 方 法 》; 所 有 的 室内 实验 均 完成 于 中 国 农业 
科学 院 水 资源 安全 高 效 利用 重点 开放 实验 室 。 
1.2.3 ”数据 处 理 

使 用 SPSS 软件 计算 统计 参数 .进行 正 态 性 检验 
和 相关 分 析 ， 采 用 地 统计 辅助 软件 GS+9.0 进行 半 方 
差 画 数 模 型 的 拟 合 与 分 维 数 的 计算 , 利用 SURFER 
软件 进行 Kriging 插值 ， 绘 制 等 值 线 图 。 计 算 原 理 可 
参见 文献 [1,25]。 


2 结果 与 分 析 
2.1 土壤 特性 及 有 效 锅 的 均值 与 变异 系数 分 析 
土壤 有 效 锅 、 含 水 率 、pH 和 有 机 质 的 均值 和 变 
异 系数 的 测定 和 计算 结果 见 表 1， 通 过 方差 分 析 书 
检验 ， 除 剖面 A 有 效 锅 含 量 的 忆 值 为 0.04， 其 余 均 
为 0, 4 项 参数 的 各 剖面 PP 值 均 小 于 0.05， 表 明 此 次 
污水 灌 没 对 土壤 有 效 锅 、 含 水 率 、pH 和 有 机 质 的 分 
布 影响 显著 。 
污水 灌溉 后 , A、B、C、D 剖面 有 效 锅 均 值 较 灌 水 
前 分 别 减 少 0.06 mgkg  、0.11 mg:kg '、0.14 mgkg 
和 0.23 mgkg :( 表 1)， 说 明 灌 水 后 一 段 时 间 内 ,土壤 


表 1 污水 灌溉 后 不 同 剖 面 土壤 特性 及 有 效 锅 的 统计 参数 表 


Table 1 Changes of soil characteristics and available cadmium contents after one sewage irrigation in different soil sections of 
sewage irrigated field 
合 水 率 pH 有 机 质 有 效 镶 
剖面 时 间 Moisture content 二 matter Available cadmium 
Section Time 平均 值 ”变异 系数 。 平均 值 ”变异 系数 其 变异 系数 平均 值 变异 系数 
Average (%) CV (%) Average CV (%) (pke) CV (%) Average (mg'kg ') CV (%) 


A 灌水 前 Before irrigation 23.18+1.40aA 6.02 8.74+0.05bB 0.59 18.7+2.60cC 3.90 0.56+0.26aA 45.53 
灌水 后 After irrigation 18.29+0.28bA 1.53 8.67+0.29aA 0.61 25.8+2.96aA 11.47 0.50+0.18aA 36.14 
B ”灌水 前 Before irrigation 12.17+3.74cdC 30.73 8.83+0.12aA 1.33 17.5+3.98dD 22.74 0.29+0.12cC 40.84 
灌水 后 After irrigation 16.37+0.44cB 2.72 8.61+0.05cC -1.13 23.2+3.98bB 7.16 0.18+0.08dC 45.16 
C 灌水 前 Before irrigation 13.19+0.94bB 7.12 8.75+0.07bB ”0.75 22.9+2.81aA 2.27 0.44+0.19bB 42.27 
灌水 后 After irrigation 18.40+0.48abA 2.63 8.69+0.10aA 0.68 26.4+3.54aA 3.41 0.30+0.11bB 36.15 
D 灌水 前 Before irrigation 12.45+0.82cBC 6.56 8.74+0.26bB 1.78 21.2+3.13bB 4.76 0.45+0.18bB 39.29 
灌水 后 After irrigation 18.50+0.92aA 4.96 8.64+0.60bB 0.59 23.4+2.92bB 2.48 0.22+0.06cC 28.31 
A、B、C 和 DD 为 研究 田 块 由 北 至 南 依次 布设 的 剖面 。 表 中 同一 时 间 内 同 列 不 同 小 写字 母 表 示 统 计 检 验 不 同 剖 面 在 P<5% 水 平 差异 显著 ; 同 
一 时 间 内 同 列 不 同 大 写字 母 表 示 统 计 检 验 不 同 剖 面 在 P<1% 水 平 差异 极 显著 。A, B, C and D refer to the sections cross the study plot from north to 


south, respectively. Different small letters in the same time and same column represent significant difference among different sections at P < 5% level. 
Different capital letters in the same time and same column represent extremely significant difference among different sections at P< 1% level. 
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垂 向 剖面 的 有 效 锅 经 过 迁移 、 植 物 吸 收 或 离子 状 
态 转 化 等 作用 而 有 所 减少 。 灌 水 前 各 剖面 有 效 镶 
变异 系数 较为 相近 ， 呈 中 等 离散 程度 ， 变 化 范围 
为 39.29%~45.53%。 灌水 后 有 效 锅 仍 呈现 出 中 等 离 
散 ， 除 剖面 B 离散 程度 增加 4.32%， 其 余 剖 面 分 别 减 
少 9.39%、6.12% 和 10.98%， 范 围 为 28.31%~ 
45.16%。 由 此 表明 ,污水 灌溉 后 土壤 中 垂 向 剖面 有 
效 锅 含量 变异 趋 于 减 小 ,有效 饮 经 过 迁移 、 转 化 等 
物理 作用 , 在 垂 向 方向 上 进行 了 重新 分 布 , 均一 化 
程度 较 明 显 。 

A、B、C、D 剖面 的 含水 率 均 值 较 灌 水 前 分 别 
变化 了 -4.89%、4.20%、5.21% 和 6.05%， 除 剖面 A 
外 ， 其 余 剖 面 土壤 合 水 率 有 所 增加 ; 位 于 地 块 北面 
的 剖面 A,， 可 能 受到 灌水 方式 (由 南 至 北 漫灌 ) 的 影 
响 ， 水 分 渗透 深度 较 南 面 土 层 浅 , 土壤 水 分 蒸发 量 
大 。 灌 水 前 剖面 B 呈 中 等 离散 程度 ， 其 余 剖 面 变 异 
系数 较为 相近 ， 变 化 范围 为 6.02%~7.12%。 灌 水 后 合 
水 率 呈 现 弱 离散 ,剖面 B 离散 程度 减少 28.01%， 其 
余 剖 面 分 别 减少 4.49%、4.49% 和 1.60%， 范围 为 
1.53%~4.96%。 由 此 表明 , 污水 灌溉 后 土壤 垂 向 剖面 
合 水 率 的 变异 减弱 。 

A、B、C、 DD 剖面 的 pH 均值 较 灌 水 前 有 所 减 小 ， 
灌水 活动 会 引起 土壤 pH 微量 变化 。 灌 水 前 后 各 剖面 
pH 均 呈 现 弱 离 散 程度 ， 除 剖面 A 离散 程度 增加 
0.02%， 其 余 剖 面 分 别 减少 0.20%、0.07% 和 1.19%。 


由 此 表明 ,污水 灌溉 后 土壤 中 垂 向 剖面 pH 的 变异 

A、B、C、D 剖面 的 有 机 质 含量 均值 较 灌 水 前 
分 别 增加 7.1 gkg '、5.7 g kg 3.5 gkg 和 2.2 gkg!; 
除 剖 面 C 的 离散 程度 增加 1.14% 外 ， 其 余 分 别 减 少 
2.43%、5.58% 和 2.28%， 呈 现 中 等 离散 程度 。 结 
显示 , 污水 灌溉 后 有 机 质 含量 有 所 增加 ， 但 分 布 的 
变异 程度 减弱 ， 表 明 此 次 污水 灌溉 后 土壤 中 垂 向 剖 
面 有 机 质 受 土壤 团 粒 结构 增多 的 影响 而 离散 程度 有 
所 减弱 。 

由 传统 统计 学 对 4 项 因子 的 分 析 数 据 可 以 看 出 ， 
污水 灌溉 后 土壤 有 效 锅 、 含 水 率 、pH 和 有 机 质 的 离 
散 程 度 有 减弱 趋势 。 

2.2 ”土壤 有 效 锅 空 间 变 异 特征 分 析 

根据 地 统计 学 的 原理 ,利用 块 金 值 (C0) 与 块 基 
比 [Cw(Cot+O)] 可 以 分 别 表征 区 域 化 变量 的 由 随机 
因素 引起 的 变异 程度 和 变量 的 空间 相关 关系 的 强 
弱 。GS+9.0 地 统计 分 析 软 件 的 半 方 差 画 数 模拟 及 
参数 值 见 表 2。 灌 水 前 后 各 剖面 的 块 金 值 差 异 较 大 ， 
灌水 后 A 剖面 的 增幅 明显 ,剖面 B、C、D 分 别 下 
降 98.6%、52.7% 和 86.2%， 变 化 范围 为 1.0x10 ~ 
7.8x10-， 说 明 随 机 性 在 剖面 A 的 空间 变异 中 占 主 
导 地 位 ， 而 污水 灌溉 后 垂 向 剖面 总 体 表 现 出 受 施 
肥 、 耕 作 措 施 、 种 植 制度 等 人 为 活动 的 影响 所 引发 
的 变异 较 小 。 


表 2 污水 灌溉 前 后 不 同 剖 面 土壤 有 效 锅 半 方差 参数 
Table 2 了 Parameters of semi-variance of soil available cadmium content in different soil sections of the sewage irrigated field before 
and after sewage irrigation 


剖面 时 间 理论 模型 块 金 值 基 台 值 块 基 比 变 程 
Section Time Theoretical model Co Co+C Co Co+C) Range (m) 
A 灌水 前 Before irrigation 球状 Spherical 1.90x10-3 6.68x10 一 0.028 0.513 

灌水 后 After irrigation 球状 Spherical 7.80x1073 4.21x10- 0.185 1.308 
B 灌水 前 Before irrigation 高 斯 Gauss 7.00x1073 2.20x10™ 0.318 1.857 
灌水 后 After irrigation 球状 Spherical 1.00x107 6.52x10 一 0.002 0.393 
C 灌水 前 Before irrigation 高 斯 Gauss 7.27x103 3.57x10 一 0.203 0.395 
灌水 后 After irrigation 高 斯 Gauss 3.44x1073 1.45x10™ 0.238 0.443 
D 灌水 前 Before irrigation 高 斯 Gauss 7.25x1073 3.62x107™ 0.200 0.419 
灌水 后 After irrigation 虽 数 Exponential 1.00x10-3 4.05x10™ 0.025 0.783 


灌水 后 Cw(Co+O) 比 值 为 0.002~0.238, 与 灌水 
前 相 比 ， 灌 水 后 剖面 A 和 C 有 所 增 大 ， 空 间 变 异 
要 由 随机 因素 引起 的 ; 剖面 B 和 D 均 大 幅 减 小 ， 随 
机 作用 减弱 。 整 体 表现 出 ,污水 灌溉 后 有 效 锅 在 垂 
向 剖面 由 随机 因素 引起 的 空间 变异 较 小 ， 呈 现 强 空 
间 相 关 性 。 

变 程 (cange) 值 呈现 出 灌水 后 土壤 有 效 锅 自 相关 


范围 为 0.393~1.308 mi; 剖面 B 的 自 相 关 范 围 变 小 ,， 
且 小 于 自 相关 范围 变 大 的 剖面 A、C 和 D， 变 程 值 (让 
表示 若 两 观测 值 之 间 的 距离 小 于 该 值 ， 则 它们 之 间 
存在 一 定 空间 相关 性 , 说 明 剖 面 B 有 效 锅 的 相关 范 
围 缩小 ; 试验 田 块 在 垂 向 剖面 上 ,土壤 有 效 锅 总 体 
呈现 空间 相关 范围 变 大 的 趋势 ， 表 现 出 一 种 小 尺度 
范围 内 的 弱 变 异 。 
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2.3 ”土壤 有 效 锅 分 维 数 分 析 

分 维 数 能 够 表示 变量 在 无 标 度 区 间 内 的 空间 
连续 程度 Pl。 图 1 为 灌水 后 各 剖面 有 效 锅 的 分 维 
数 图 ， 由 图 1 可 以 看 出 ,灌水 后 各 剖面 的 分 维 数 均 
处 于 1~2 间 ， 表示 着 土壤 中 饥 的 有 效 态 参数 值 在 
连续 面 上 处 于 严格 的 空间 依赖 和 完全 的 空间 独立 
之 间 , 土壤 结构 具有 分 形 特征 ， 两 者 是 具有 内 
在 联系 的 。 
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图 1 


试验 中 灌水 前 ,剖面 A、B、C 和 D 的 分 维 数 分 
别 为 1.883、1.776、1.906 和 1.887， 灌 水 后 的 土壤 有 
效 态 包 的 半 方 差 画 数 [y(h)]- 变 程 (h) 关 系 县 有 分 形 规 
律 , 如 图 1 所 示 。 剖 面 B 的 有 效 饥 分 维 数 最 高 ， 为 
1.997; 其 余 剖 面 的 分 维 数 均 较 灌水 前 降低 ， 分 别 为 
剖面 A1.747、 剖 面 C 1.875 和 剖面 D 1.840， 总 体 上 
表明 有 效 锅 在 垂 向 分 布 上 的 空间 变异 由 结构 性 因素 
起 决定 作用 。 
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污水 灌流 后 剖面 A、B、C 和 D 的 土壤 有 效 锅 的 分 维 数 


Fig.1 Fractal dimensions of soil available cadmium in sections A, B, C and D in sewage irrigated field after irrigation 
六 半 方 差 画 数 ; h: 变 程 。x: semivariance; h: separation distance. 


土壤 有 效 锅 垂 向 分 布 分 析 

等 值 线 图 可 以 直观 和 准确 地 描述 试验 区 不 同 剖 
面 土壤 有 效 态 锅 的 垂 向 分 布 情况 。 试 验 田 块 为 北 东 
高 ， 南 西 低 ， 粉 砂 质 黏土 ; 灌水 方向 为 由 南 向 北 。 试 
验 田 块 依据 试验 数据 ， 绘 制 的 等 值 线 图 如 图 2 所 示 。 
从 图 中 可 以 看 出 ,灌水 后 , 剖面 A.B 和 C 土壤 有 效 
态 锅 大 体 呈 现 东 西向 的 条 带 状 分 布 ,剖面 D 呈现 为 
块 状 分 布 。 剖 面 上 土壤 有 效 锅 均 在 表层 0~15 cm 之 
间 含 量 较 高 ,表现 出 一 定 的 表 聚 现象 ; 较 低 含量 值 
出 现在 25~35 cm 的 范围 内 。 由 剖面 A 至 D， 有 效 饥 
含量 等 值 线 分 布 由 密集 到 稀疏 ， 变 化 幅度 逐渐 变 小 ， 
梯度 减弱 ， 方 向 有 所 改变 。 漫 灌水 流 由 A 至 D 过程 
中 , 因 地 势 和 土壤 颗粒 大 小 等 因素 的 影响 ， 表 现 出 


2.4 


有 效 锅 含 量 在 垂 向 剖面 上 由 复杂 至 简单 的 分 布 差异 ， 
说 明 有 效 锅 合 量 的 垂 向 空间 分 布 主要 由 结构 性 因素 
控制 。 

各 剖面 土壤 有 效 饥 含 量 灌水 后 均 有 所 减少 ,等 
值 线 梯度 减弱 ， 趋 向 于 块 状 分 布 ， 表 明 灌 水 之 后 ， 
土壤 有 效 锅 发 生 迁 移 和 转换 ,同时 水 中 的 锅 也 会 受 
到 土壤 颗粒 的 吸附 ， 形 成 一 定 累积 效应 。 土 壤 有 效 
锅 的 分 布 受 到 灌水 影响 ,表现 出 剖面 上 的 均 质 分 布 ， 


变异 性 减弱 。 
2.5 ”土壤 特性 及 有 效 锅 的 相关 分 析 


统计 学 中 ,相关 系数 可 以 用 来 表征 变量 间 是 否 
存在 某 种 相关 关系 及 相关 的 强 弱 程度 。 表 3 为 土 
壤 含 水 率 ` 有机质 和 pH 与 土壤 有 效 锅 含量 相关 性 结 
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2 污染 灌溉 前 后 剖面 A、B、C、D 土壤 有 效 态 锅 的 等 值 图 
Fig. 2 Contour maps of soil available cadmium in section A, B, C and D in sewage irrigated field before and after sewage irrigation 


表 3 污水 灌流 不 同 剖 面 土壤 有 效 锅 与 影响 参数 的 相关 关系 
Table 3 “Correlation coefficients between soil available cadmium and influencing factors in different sections of sewage irrigated 
field before and after sewage irrigation 


影响 参数 Influencing factor 时 间 Time A B C D 
土壤 合 水 率 灌水 前 Before irrigation 0.876™ 0.273™ -0.424™ -0.833™ 
Soil moisture 灌水 后 After irrigation 0.405™ 0.511™ 0.752™ 0.819” 
土壤 有 机 质 灌水 前 Before irrigation 0.907™ 0.593™ 0.944” 0.907™ 
Soil organic matter 灌水 后 After irrigation 0.533” 0.782™ 0.937™ 0.573™ 
pH 灌水 前 Before irrigation -0.745” -0.607™ -0.669™ -0.907™ 
灌水 后 After irrigation 0.024 -0.572™ -0.705™ -0.640” 


* 表 示 Sig.( 双 侧 )<0.05, 呈 显 著 相 关 ; ** 表 示 Sig.( 双 侧 )<0.01, 呈 极 显著 相关 。* means significant correlation at Sig. < 0.05 level; ** means 
significant correlation at Sig. < 0.01 level. 
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果 。 从 表 3 可 以 看 出 ， 有 效 锅 与 有 机 质 呈 强 相 关 ， 与 
土壤 含水 率 呈 中 等 程度 相关 , 与 pH 呈现 强 负 相关 。 
土壤 有 机 质 的 含量 越 高 ， 土壤 有 效 锅 含 量 越 大 ; pH 
越 高 ， 土壤 有 效 锅 含 量 越 小 。 灌 水 之 后 , 土壤 含水 率 
与 有 效 锅 的 相关 程度 变化 较 大 ， 灌 水 会 削弱 有 效 锅 
与 有 机 质 和 pH 的 相关 程度 ,说 明 土 壤 中 有 效 锅 合 
量 受 到 土壤 有 机 质 、 离 子 交换 作用 、 灌 溉 水 质 以 及 
土壤 环境 中 pH 等 因素 的 影响 。 
3 讨论 与 结论 

本 研究 采用 地 质 统计 学 与 传统 统计 学 相 结合 
的 方式 ， 对 灌区 小 尺度 田 块 单 次 污水 灌溉 垂 向 剖面 
土壤 有 效 锅 含 量 的 空间 分 布 和 变异 进行 分 析 ， 结 
表明 : 

污水 灌溉 后 ,土壤 有 效 馈 合 量 较 灌水 前 有 不 同 
程度 减少 ， 空 间 分 布 呈 中 等 离散 程度 ,土壤 含水 率 
和 土壤 pH 呈现 弱 离 散 ， 灌 水 会 减弱 土壤 有 效 锅 、 含 
水 率 、pH 和 有 机 质 的 空间 离散 程度 。 考 虑 到 灌溉 水 
中 含有 锅 ， 灌 水 后 田 块 各 剖面 的 土壤 有 效 锅 含 量 
未 如 预期 一 样 有 所 增加 ， 说 明 在 溶质 运 移 的 作用 下 
土壤 中 重金 属 向 下 迁移 C)， 淋 溶 和 迁移 的 含量 大 于 
水 体 向 土壤 输入 的 含量 。 

灌水 前 土壤 有 效 锅 的 空间 变异 受 随机 因素 的 影 
响 较 大 ， 这 与 文献 [11] 中 土壤 锅 的 结论 较为 相近 ,但 
污水 灌 涉 后 , 土壤 有 效 久 分 维 数 变 小 ，Co/(Co+C) 剖 
面 间 呈 波 动 变化 ， 垂 向 剖面 总 体 受 随机 性 因素 影响 
较 小 ， 空 间 变 异 程度 弱 ， 自 相关 性 强烈 ， 剖 面 间 的 
空间 变异 性 与 田 块 总 体 变化 不 一 致 ， 其 波动 变化 的 
情况 受 人 为 、 土 壤 结 构 和 地 域 环境 的 相互 作用 ， 影 
响 原因 需要 进一步 试验 和 探讨 。 分 维 数 与 地 统计 方 
法 虽然 在 理论 上 具有 一 定 的 数学 关系 和 意义 ， 但 在 
本 次 研究 中 两 种 方法 的 试验 结果 并 不 完全 统一 ， 灌 
水 后 , 剖面 B 和 DD 的 土壤 有 效 锅 分 形 和 变异 与 理论 
推断 出 现 差 异 ， 这 可 能 是 由 于 土壤 并 非 是 具有 理想 
分 形 特 征 的 介质 2 土壤 有 效 锅 基 于 土壤 介质 的 分 
形 研究 也 有 待 进一步 验证 。 但 分 维 数 理论 的 分 析 结 
果 仍 可 以 说 明 在 小 尺度 田 块 中 , 土壤 有 效 锅 也 适宜 
采用 地 统计 模型 进行 空间 分 析 和 与 估 值 "1。 

等 值 线 分 布 显示 ,各 剖面 土壤 有 效 锅 垂 向 空间 
分 布 呈 现 条 带 状 ， 受 地 势 和 灌水 方向 的 影响 ,各 剖 
面 有 不 同 分 布 情况 ,但 表土 层 (0~15 cm) 合 量 均 大 于 
心 土 层 (30~40 cm)， 可 见 锅 元 素 主 要 富 集 在 土壤 表 
层 ,这 与 大 多 数 研究 者 的 结论 较为 一 致 "土壤 有 
效 锅 纵 向 不 均匀 性 可 能 受到 土壤 理化 性 质 和 植物 根 
系 吸 收 的 共同 作用 。 


另外 , 土壤 有 效 锅 合 量 与 土壤 有 机 质 含量 呈正 
相关 , 与 土壤 含水 率 呈 中 等 程度 相关 , 与 pH 呈现 强 
负 相 关 ， 说 明 土 壤 pH、 有 机 质 含 量 、 黏 粒 的 含量 可 
以 影响 土壤 中 重金 属 的 含量 后 )， 同 时 土壤 对 铀 的 吸 
附 可 能 主要 取决 于 土壤 pH 水 平和 碳酸 盐 含 量 广 ]; 
污水 灌溉 对 各 因子 的 影响 显著 ， 它 会 削弱 有 效 饥 与 
影响 因子 间 的 相关 关系 ， 从 而 影响 土壤 有 效 锅 的 
空间 变异 。 
研究 发 现 4 块 剖面 土壤 有 效 锅 分 布 和 变异 基本 
可 以 反映 出 田 块 整体 有 效 锅 的 情况 ; 单 次 污水 灌溉 
扰动 后 ， 有 效 锅 在 水 和 土壤 中 相互 置换 ,土壤 中 有 
效 锅 发 生 迁 移 和 再 分 布 ， 进 而 在 垂 向 剖面 上 形成 空 
间 变 异 ; 灌水 作为 连通 土壤 系统 和 植物 系统 的 关键 
因素 ， 可 以 认为 植物 吸收 作用 对 于 土壤 有 效 锅 的 分 
布 和 变异 是 不 可 忽视 的 。 若 能 结合 修复 植物 或 作物 
的 样本 数据 ， 分 析 该 田 块 土壤 有 效 锅 在 土壤 -植物 
系统 中 吸附 和 迁移 的 初步 关系 ， 对 于 土壤 重金 属 污 
染 植 物 修复 方面 的 研究 县 有 一 定 意义 。 
由 于 土壤 重金 属 垂 向 剖面 分 布 和 变异 的 研究 方 
式 受 采样 深度 和 广度 的 限制 ， 此 次 研究 结论 不 能 较 
好 地 适用 于 大 尺度 范围 的 试验 分 析 。 另 外 ,研究 中 
的 灌水 动态 仅 为 单 次 , 而 由 污水 灌溉 引起 的 土壤 重 
金属 污染 是 一 个 长 期 渐变 的 过 程 ， 需 要 在 水 质 、 灌 
水 时 间 和 频次 方面 做 进一步 的 探讨 。 
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